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We show the existence of a system size coherence resonance effect for an ensemble of globally coupled 
excitable systems. Namely, we demonstrate numerically that the regularity in the signal emitted by an ensemble 
of globally coupled Fitzhugh-Nagumo systems, under excitation by independent noise sources, is optimal for a 
particular value of the number of coupled systems. This resonance is shown through several different dynamical 
measures: the time correlation function, correlation time and jitter. 
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Noise induced resonance is a topic that has attracted a lot of attention in the last years. In particular, it has been unambiguously 
shown that the response of some systems to an external perturbation can be enhanced by the presence of noise (stochastic 
resonance [[!]-§]]). A different effect is that of coherence resonance [||| in which an excitable system shows a maximum degree 
of regularity in the emitted signal in the presence of the right amount of fluctuations (or the related one of stochastic resonance 
£f~) • without the need of an external forcing Coherence resonance has also been studied in dynamical systems close to the 
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neuronal models such as the FitzHugh-Nagumo [|11J,|12|] and Hodgkin-Huxley J13Q models. It has been observed experimentally 
in electronic circuits, either excitable Jl4l , |l5[ ] or chaotic [ |l6||l^ ], and in lasers operating in an excitable regime [|l8|]. 

In an important recent paper fll^], Pikovsky et al. have shown that when one considers an ensemble of coupled bistable 
lf~) ' systems subjected to an external periodic forcing (and in the presence of a constant amount of noise), it turns out that an optimal 
response is obtained for an appropriate value of the number N of coupled systems. In other words, that there is a resonance 
with respect to the number of coupled elements, rather to the usual one that involves the noise level. The authors speculate that 
this system size resonance might be relevant to neuronal dynamics, in which the neuronal connections or the coupling strengths 
between neurons can be tuned in order to achieve maximum sensitivity to external signals. 

In this paper, we extend the previous result by considering an ensemble of globally coupled excitable systems, each one under 
the influence of its own noise with a fixed intensity, but without an external forcing. We show that there is a coherence resonance 
effect as a function of the number N of coupled systems. More specifically, we show that the excitable systems pulse on average 
with a regularity which is optimal for a specific value of N. 



Motivated by the biological applications suggested in [19] we consider an ensemble of i = 1, . . . , N coupled FitzHugh- 



Nagumo systems, each one described by the activator, Xi, and inhibitor, y.;, variables. The FitzHugh-Nagumo model provides 
the simplest representation of firing dynamics and has been widely used as a prototypic model for spiking neurons as well as for 
cardiac cells po| , ^l| ]. The dynamical FitzHugh-Nagumo equations modified to account for the global coupling are as follows: 
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ex\ = Xi - -x 3 - yi + — VVxj ■- Xi) (1) 

: i=i 

yi = x i + a + D£ i (t) (2) 

where independent noises of intensity D have been added to the slow variables y$ as in ref. [pp. The are white noises with 
Gaussian distribution of zero mean and correlations (£,i(t)£,j{t')) = 5ijS(t — t'). The difference in the time scales of Xi and yi 
is measured by e, a small number. The systems are globally coupled, as indicated by the last term of Eq. ([!]), where K is the 
coupling strength. 

We consider the excitable regime, a > 1. In the absence of coupling, K = 0, each FitzHugh-Nagumo system emits pulses. 
The pulses are trajectories that exit the basin of attraction of the stable fixed point and are triggered by the influence of the noise 
term. For zero coupling, coherence resonance exists when the regularity of the time between pulses is optimal for a certain value 
of the noise intensity D [||]. This behavior is the same for all the systems although the response of each one is uncorrelated with 
any other. 
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Let us now study the collective response of the coupled system and compare it with the individual responses. For the collective 
response, we introduce the average values of the activator and inhibitor variables as 



N N 

By following the approach by Desai and Zwanzig [[22j| (see also reference JI^]) it is possible to reach an approximate effective 
equation for these average values of the form: 

eX = F(X) - Y (4) 

Y = X + a + -=£(i) (5) 
v N 



where is a white noise source. Although the exact form of the function F(X) and the analysis of the approximations 
assumed in the derivation will be presented elsewhere, we need only to remark here that in the (exact) equation for Y(t) the 
noise intensity appears rescaled as D/V~N. Therefore, this approximation suggests that the optimal effective noise intensity 
for the appearance of coherence resonance can be achieved by varying the number of coupled elements N, as in the case of 
stochastic resonance for the bistable system considered in [|l9[]. To go beyond this approximation, we numerically integrate the 
equations of motion (1) and (2). 

Figure [j] (left panel) shows the time trace for the variable X(t) while |l] (right panel) shows the time trace for the variable Xi (t) 
of one of the elements chosen randomly, for three different values of the number of coupled elements (see the caption of the 
figure for details of the parameters). Notice that for N = 160 the regularity of the emitted pulses is better than that corresponding 
to larger or smaller values of N. This is a clear signature of coherence resonance. Moreover, it can be seen that the regularity in 
the averaged variable X(t) is better than in one of the individual elements, showing that the coupling allows for a smoothness of 
the trace. It is worth noting that the peaks in the collective variable X(t) and Xi(t) are very well synchronized in time indicating 
that the individual systems are pulsing synchronously in time. In Figure || (left panel) we plot the time trace for the slow variable 
Y(t), as well as a time trace for a single yi(t) (right panel). At variance with the fast variable X, it turns out that the averaged 
Y(t) shows a very nice regular behavior for an intermediate number of elements while the individual traces yi(t) do not. 

We have computed two indicators commonly used to quantify this effect [gj], First, we have computed the time correlation 
function Cx if) of the averaged X variable, defined as 

gvw= <"$y» ««)=x (t , -<*(o> (« 

and similarly for the correlation function Cy (t) for the averaged Y variable. Here the averages ( ) are with respect to the time 
t' , after a small transient has been neglected. Figure || shows this correlation function for both the X and Y variables. It can be 
seen that the correlations extend further in time for an intermediate value, neither very large nor very small, of the number of 
coupled systems N. To obtain a quantitative indicator of this effect, we define the characteristic correlation times tx and Ty for 
each variable as 

/•OO 

tx.y = / \C x ,y(t)\dt (7) 
Jo 

In practice, the upper limit of the integral is replaced by a value t max such that the correlation function can be considered as 
decayed to its asymptotic value Cx.y = (t rnax — 50 for the data shown in figure fy. We have plotted these two correlation 
times in the left panel of figure Q Both times reach a maximum at approximately the same value N w 160, indicating that, for 
the set of parameters chosen, the maximum extent of the time correlation occurs for this number of coupled excitable systems. 

Another common indicator for the regularity of the emitted pulses can be obtained by the jitter of the time between pulses [^]. 
A pulse in the X(t) variable is defined when X(t) exceeds a certain threshold value Xq (taken arbitrarily as Xq = 0.3, although 
other values yield similar results). The jitter Rx is defined as the root mean square of the time Tx between two consecutive 
pulses normalized to its mean value: 

Rx - w (8) 

and an equivalent definition for the jitter Ry of the Y variable. The smaller the value of Rx.y, the larger the regularity of the 
pulses (a value of Rx.y = indicates a perfectly periodic signal). It is shown in the right panel of figure^] that indeed the jitter 
in both variables have a well defined minimum at a value of N w 80, again showing the existence of the system size resonance. 
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When comparing with the results of the correlation time, it is not uncommon that the two indicators (the correlation time r and 
the jitter R) have their optimal values at different values of the system parameters [||Jl^]. 

In summary, we have shown that an ensemble of globally coupled FitzHugh-Nagumo excitable systems subjected to inde- 
pendent noises pulse on average with a regularity that is maximum for a given value of the number N of coupled systems. An 
approximate calculation indicates that the collective variable Y(t) is subjected to a noise of effective intensity D/y/N. There- 
fore, even in the presence of a large amount of noise (D large), it is possible to couple the right number of systems in order to 
optimize the periodicity of the emitted pulses. Since the FitzHugh-Nagumo system has been used previously to model some 
biological systems, we believe that our results can be relevant when analyzing the collective response of such systems in a noisy 
environment and can help to explain the observed size of some ensembles of excitable cells in living organisms. 
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FIG. 1. Time series for the averaged variable X{t) (left panel), and for the individual variable xi(t) (right panel) of the set of coupled 
FitzHugh-Nagumo systems, as obtained from a numerical integration of Eqs.(|l|-^|), for different values of the number of coupled elements: 
N — 1 (top), N — 160 (middle) and N = 1000 (bottom). Observe that the largest regularity is obtained for the intermediate value of N. 
The equations have been integrated numerically using a stochastic Runge-Kutta method (known as the Heun method [23]) with a time step 
h = 10 -4 and setting the following parameters: a = 1.1, e = 0.01, K = 2, D = 0.7. 
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FIG. 2. Time series for the averaged variable Y(t) (left panel), and the individual variable yi(t) (right panel) of the set of 
FitzHugh-Nagumo systems, Eqs.(|l|-^|). Similarly as in Figure [IJ observe that again the largest regularity for the averaged Y variable is 
obtained for the intermediate value of N. In this case, however, there is no obvious increase in the regularity of the yi individual variables. 
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FIG. 3. Correlation functions Cx(t) and CV(t) of the averaged variables X(t) and Y(t), respectively, for the cases of TV = 1 (dotted 
line), TV = 160 (solid line) and TV = 1000 (dashed line). Notice that, in agreement with the qualitative results derived from figures [l] and ^, 
the slower decay of the correlations corresponds to the intermediate values of the system size TV. Same parameters as in figure hi 




FIG. 4. Panel (a) plots the correlation times tx and ry as obtained by integration of the absolute value of the respective correlation 
functions. Clear maxima (maximum extent of the correlations) can be observed around TV = 160. Panel (b) plots the jitter of the time between 
consecutive pulses of the collective variables X(t) (stars) and Y(t) (triangles). Clear minima (optimal regularity in the emitted pulses) can be 
observed around TV = 80 in both cases. 
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